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Summary 

The  resi  arch  performed  under  the  contract,  during  the  period  26  September  1991 
through  25  Sep  ember  1992,  can  be  divided  into  two  main  topics;  using  the  moment  of 
regional  Raylei  gh  waves  and  a  regional  magnitude  based  on  the  maximum  amplitude  of 
the  shear  wave  train  as  a  seismic  discriminant  and  source  retrieval  from  broadband 
regional  seismograms. 

In  section  1,  we  propose  using  the  ratio  of  Ml  (local  magnitude)  to  M()  (scalar 
seismic  moment)  as  a  regional  discriminant  analagous  to  the  teleseismic  discriminant 
using  the  m^tM ,  ratio.  We  applied  this  criterion  to  a  data  set  of  299  earthquakes  and  178 
explosions  and  found  that  this  ratio  appears  to  be  diagnostic  of  source  type.  For  a  given 
M„,  the  Ml  of  an  explosion  is  more  than  0.5  magnitude  units  larger  than  that  of  an 
earthquake.  This  separation  of  populations  with  respect  to  source  type  can  be  attributed 
to  the  fact  that  Ml  is  a  short-period  (I  11/)  energy  measurement,  whereas  seismic  moment 
is  determined  Iront  long-period  body  wave  phases  (period  >  4  seconds)  and  surface 
waves  ( 10  to  40  seconds).  Using  regional  stations  with  sources  200  to  600  km  away,  the 
effective  threshold  for  magnitude  measurements  for  this  discriminant  is  found  to  be  Ml  = 
3.1  for  earthquakes  and  Ml  3.6  for  explosions.  This  method  does  require  the 
determination  of  regional  crustal  models  and  path  calibrations  from  master  events  or  by 
other  means. 


In  section  2,  we  develop  and  test  a  method  of  relocation  and  source 
characterization  of  small  earthquakes  using  one  modern  regional  station.  First,  we  model 
teleseismic  body-waves  of  two  events,  which  are  used  as  masters.  Short  period  depth 
phases,  pP  and  ;P,  are  used  to  establish  the  epiccntral  depth,  and  the  events  are  relocated 
using  calibrated  stations  and  a  mantle  model  derived  for  this  region,  TIP.  The  events 
moved  upward  by  12  and  31  kms,  respectively.  The  regional  waveforms  recorded  at  the 
IRIS  station  GAR  from  the  best  determined  source  are  forward  modeled  to  establish  a 
local  crustal  mi  del.  A  four-layered  model  with  a  thickness  of  65  km  proves  effective. 
Synthetics  from  this  crustal  model  can  then  be  compared  to  the  data  from  other  events 
where  the  depths  arc  poorly  known,  i.c.,  PDF  depths  of  33  km.  Forward  modeling  of  the 
short  period  phases  allows  better  estimates  of  depth  and  makes  relocation  possible.  Next, 
we  perform  a  long  period  inversion  of  whole  and  partial  waveform  data  to  obtain  source 
mechanism  and  moment,  and  repeat  the  procedure  if  necessary.  This  technique  is  applied 
to  events  in  the  tectonically  active  Pamir-Hindu  Kush  region  to  test  its  usefulness.  Five 
crustal  events  sampling  various  azimuths  are  presented  as  examples  of  relocation  and 
determination  of  source  mechanisms  from  small  events.  The  smallest  event  has  a 


moment  of  5.1  x  10“  dyne-cm.  For  the  seven  events  we  studied  the  average  depth 


correction  to  PI  E  is  19  km,  and  the  location  correction  is  1 1  km  on  average.  The  method 
can  he  used  to  identify  earthquakes  and  thus  it  lowers  the  threshold  below  present  mt,:Ms 
discrimination,  since  these  events  are  too  small  to  be  seen  teleseismically.  Secondly, 
these  events  cu  i  be  used  as  "masters”  in  the  calibration  of  other  systems,  case-based 


event  character!  nations,  etc. 
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MliMq  as  a  Regional  Seismic  Discriminant 
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Mi’.Mq  as  a  Regional  Seismic  Discriminant 

Bradley  B.  Woods,  Sharon  Kedar  and  Donald  V.  Helmberger 
Seismological  Laboratory  852-21, 

California  Institute  of  Technology 

The  mg; Ms  ratio  determined  by  teleseismic  observations  has  proven  to  be  an  effective 

discriminant,  for  explosive  sources  tend  to  be  significantly  richer  in  short-period  energy  than 
earthquakes.  Unfortunately,  this  method  is  limited  by  the  detection  threshold  of  teleseismic 
surface  waves.  However,  recent  advances  in  instrumentation  allowing  low  ampliiude  surface 
wave  measurements  coupled  with  new  analytical  techniques  make  it  feasible  to  use  regional 
wa  veform  data  to  determine  the  long-period  source  excitation  level  of  low  niagnii  ude  events. 
We  propose  using  the  ratio  of  A//,  (local  magnitude)  to  M0  (scalar  seismic  moment)  as  an 
ana  logouts  regional  discriminant.  We  applied  this  criterion  to  a  data  set  of  299  •  arthquakes 
and  176  explosions  and  found  that  this  ratio  appears  to  be  diagnostic  of  source  type.  For  a 
given  V/n,  the  A//,  of  an  explosion  is  more  than  0.5  magnitude  units  larger  thai  that  of  an 
earl  liqunke.  This  separation  of  populations  with  respect  to  source  type  can  hi  attributed 
to  the  fact  that  Af/,  is  a  short-period  (1  Hz)  energy  measurement,  whereas  seismic  moment 
is  determined  from  long- period  body  wave  phases  (period  >  4  seconds)  and  surface  waves 
(10  to  40  seconds).  Using  regional  stations  with  sources  200  to  600  km  away,  the  effective 
threshold  for  magnitude  measurements  for  this  discriminant  is  found  to  be  Mi,  =  3.1  for 
earthquakes  and  Ml  -  3.6  for  explosions.  This  method  does  require  the  determination  of 
regional  crustal  models  and  path  calibrations  from  master  events  or  by  other  means. 
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Introduction 


In  this  era  of  increased  concern  regarding  t,he  proliferation  of  nuclear  weapons,  the  need 
for  effective  seismic  discrimination  techniques  is  as  important  as  ever.  As  more  countries 
attain  or  near  the  technology  threshold  needed  to  develop  nuclear  weapons,  more  regions  of 
the  world  need  to  he  monitored  for  verification  purposes.  Effective  discrimination  methods 
which  can  make  u-;o  of  historical  seismic  data  and  don’t  require  elaborate,  costly  large-scale 
arrays,  are  of  com iderable  interest. 

One  of  the  most  successful  seismic  discriminants  proven  so  far  is  the  classical  comparison  of 
body  wave  magnitude  (mi,)  to  surface  wave  magnitude  (Ms)  (Basham,  1%!);  Lieberman  and 
I’oincroy.  1909:  and  Marshall,  1 1)7 1 :  Stevens  and  Day,  1985).  which  exploits  the  observation 
that  for  a  given  no,.  <  x plosions  have  a  significantly  smaller  Ms  than  do  earthquakes.  This 
observation  implies  that  explosive  sources  are  richer  in  high  frequency  energy  than  are  earth¬ 
quakes  for  a  given  long-period  energy  level  and  is  attributed,  in  part,  to  the  differences  in 
characteristic  temporal  and  spatial  source  dimensions  between  the  two  source  types.  Savino 
ft  a!  (1971),  Aki  it  al  (1974)  and  Muller  (1973)  find  that  explosions  exhibit  characteristics 
of  an  impulsive  source.  Empirical  explosion  source  models  developed  by  Haskell  (1967), 
and  M filler  and  Murphv  (1971)  yield  source  time  functions  with  rise'  times  which  are  only 
fractions  of  second  in  length.  Aki  (1967).  Bruno  (1970)  and  Marshall  (1970)  find  that,  earth¬ 
quake  source  fund  ions,  on  the  other  hand,  are  best-modeled  as  ramps  or  step  functions  with 
combined  duration  and  rise  times  greater  than  one  second.  Dreger  and  Ilelinberger  (1991a) 
found  that  even  moderate-sized  earthquakes  (Mi.  ~  5.2)  have  source  durations  greater  than 
a  second.  Stevens  and  Day  (1985)  conclude  from  numerical  modeling  experiments  that  the 
difference  between  earthquake  and  explosion  source  spectra  is  only  partially  responsible  for 
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explaining  observations.  They  also  cite  focal  mechanism,  near-source  elastic  proper¬ 

ties  and  pP  interference  effects  as  contributing  factors  to  the  separation  of  populations  with 
this  type  of  discriminant. 

The  drawback  to  the  mt,:Ms  method  is  the  threshold  at  which  teleseismic  surface  wave 
magnitudes  can  be  determined  and  the  apparent  convergence  of  populations  ai  small  magni¬ 
tudes.  Lieberman  and  Pomeroy  ( I960)  found  the  surface  wave  detection  level  for  earthquakes 
to  be  rn b  >4.3  and  for  explosions  to  be  m<,  >  4.8,  but  found  that  the  populations  converged 
below  iub  =  5.0  for  discrimination  purposes.  Evernden  tl  al  (1971)  suggest  that  with  high- 
dynatnic-range  digital  instruments,  surface  waves  from  earthquakes  with  m,,  >4.0  (which 
corresponds  roughly  to  a  m*  =5.0  explosion)  can  be  measured  out  to  6000  km  to  7000  km  - 
di  ‘ances  at  which  well-dispersed  20  sec.  surface  waves  can  be  measured.  Broadband,  high- 
dyriamic- range  seismic  stations,  such  as  those  in  the  IRIS  network  achieve  this  observational 
capability. 

Below  these  threshold  levels  the  surface  wave  signals  are  within  the  noise  lecel  and  analo¬ 
gous  discriminants  using  regional  phases  must  be  employed.  Evernden  et  al  (1971),  Lambert 
and  Alexander  (1971),  and  Peppin  and  McF.villv  (1974)  found  that,  they  con'd  distinguish 
between  source  types  on  the  basis  or  regional  Pn  amplitude  (or  »«(,(/>„))  to  regional  Airy  phase 
Havleigh  wave  amplitude  (or  magnitude)  comparisons  for  events  down  to  m,jpnj  =  3.6  to 
4.0  -  the  detection  threshold  for  explosion  generated  surface  waves. 

Two  significant  advances  in  observational  seismology  occurred  since  these  siudies:  one  is 
the  large-scale  upgrade  and  augmentation  of  seismic  networks  with  broadband,  iiigh  dynamic 
range  instruments;  the  other  is  the  increased  sophistication  in  source  parameterization  of 
earthquakes  and  explosions  in  terms  of  the  excitation  of  Green’s  functions  ar  d  the  seismic 
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moment  tensor  s»-lutions( Dziowonski  ft  til,  1981),  wiiich  more  precisely  quantify  seismic 
sources.  With  several  three  component,  broadband  stations  it.  is  feasible  to  invert  for  the 
source  function  using  regional  body  waves  (Dreger  and  Helmberger,  1990  &  1992),  surface- 
waves  (Kanamori  and  (liven.  1981;  Thioand  Kanamori,  1992;  Patton,  1988)  or  a  combinat  ion 
< >1  i  iie  t  wo  ware  t\  :>os  (Zhao  and  Helmberger,  199.1).  Mo  should  better  reflect  the  long-period 
source  characters!  ics  of  a  source  than  does  Ms,  so  that  its  use  should  improve  discriminants 
which  make  use  of  long-period  seismic  phases. 

We  re-examine  the  utility  of  short-period  vs.  long-period  energy  measures  in  the  context 
of  regional  seismic  records.  Figure  1  c' .splavs  the  study  area,  which  includes  central  and 
sout  hern  California,  western  Nevada  and  northern  Baja,  Mexico.  Waveform  data  is  presented 
in  this  paper  for  the  events  shown  in  the  figure  (stars).  This  region  is  unique  in  that  there 
is  a  large  amount  of  natural  seismicity,  as  well  as  that  of  Nevada  Test  Site  (NTS)  explosions 
and  their  aftershocks.  The  large  number  of  seismic  networks  monitoring  this  regions  has 
created  a  wealth  of  data  useful  for  discrimination  studies. 

Wit  h  Mo  servin'-  ?*  a  long-period  energy  measure,  we  employ  the  classical  local  magnitude, 
.1//..  I  Richter,  193d  &  1958)  as  a  measure  of  short-period  energy.  M /,  is  computed  from  the 
peak  horizontal  d  splacements  on  a  Wood- Anderson  torsion  instrument,  which  is  a  high 
pas-  filter  peaked  around  1.0  second.  Teles  -’smic  P-wave  phases  measured  to  determine 
body-wave  magnitudes  are  typically  of  the  same  period,  so  Mi,  too,  can  be  considered  a  fair 
measure  of  the  sho  t-period  source  spect  rum.  There  are  several  advantages  of  using  Mi.  One 
is  that  it  is  a  simple  measurement  to  make  and  is  easy  to  obtain  to  very  small  magnitudes. 
Secondly,  it  is  a  routine  source  quantification  used  by  many  seismic  networks.  We  will  make 
use  of  Mi  magnitudes  from  several  network  catalogs.  One  drawback  to  measuring  Ml  is  that 
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it  is  not  determined  from  a  particular  seismic  phase.  Normally  it  is  a  measure  of  shear-wave 
amplitude,  but  depending  on  distance  and  source  spectrum  for  extreme  cireumstances,  a 
P-wavc  or  surface  wave  may  be  the  largest  amplitude  phase  on  a  short-period  record. 

From  visual  inspection  of  broadband  recordings  of  regional  earthquakes  and  explosions  it 
is  apparent  that  these  two  source  types  show  distinctly  different  spectral  content.  Figure 
2a  compares  records  of  the  explosion  Kearsarge  (the  150  kt,  =  5.5,  .Joint  Verification 
Event  detonated  at  Pahute  Mesa)  to  two  earthquakes:  Lee  Vining  (10/24/90,  m*  =  5.0) 
and  Little  Skull  Mountain  (6/29/92,  mt  =  5.7).  The  three  events,  all  recorded  digitally 
at  Pasadena,  California  (PAS),  have  roughly  the  same  epicentral  distance  (see  Figure  1). 
For  each  event  the  broadband  displacement  records  are  displayed  (middle  trace),  as  well 
as  the  displacements  convolved  with  a  Press-Ewing  30-90  (PE)  long-period  instrument  (top 
trace)  and  with  a  Wood-Anderson  short-period  (WA)  instrument  (bottom  trace).  Peak 
amplitudes  are  to  the  right  of  each  trace.  For  Eec-sarge,  the  ratio  of  peak  short-period 
(WA  trace)  to  long-period  (PE  trace)  amplitude  is  0.325,  Q.265  and  0.625  for  the  vertical, 
radial  and  tangential  component,  respectively.  On  the  long-period  tangential  component,  no 
fundamental  Love  wave  is  apparent;  only  later  arriving,  higher  frequency,  higher  mode  waves 
(/>,,  and  possibly  imilti-pathing  crustal  waveguide  surface  waves)  are  evident.  The  absence 
of  a  hove  wave  is  also  evidence  of  an  isotropic  source,  since  explosions,  in  the  absence  of 
leclonic  strain  release  do  not  generate  long-period  Sll  waves.  For  Lee  Vining  the  ratio  of 
peak  short-period  to  long-period  amplitude  is  0.152,  0.223  <  nd  0.210  for  the  vertical,  radial 
and  tangential  components.  The  tangential  component  also  has  a  large  fundamental  mode 
Love  wave.  For  Little  Skull  Mountain  the  ratio  of  peak  short-period  to  long-period  amplitude 
is  0.0451, 0.0433  and  0.0872  for  the  vertical,  radial  and  tangential  components.  Again  a  large 
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Love  wave  is  observed.  Of  these  three  events,  the  explosion  has  the  largest  short-period  to 
long  period  am  pi  i  Hide  ratio  for  each  component. 

!■  igure  *21)  shows  an  analogous  plot  of  seismograms  recorded  at  the  digital  station  Goldstone, 
California  (CSC),  for  the  explosion  Bexar  (mi  =5.6  and  detonated  within  a  kilometer  of  the 
Kearsarge  shot  point)  and  the  two  earthquakes  of  Figure  2.  As  with  Kearsarge,  little  or  no 
long  period  Love  wave  is  evident  on  the  tangential  component;  only  shorter  period  crustal 
waveguide  surface  waves  are  visible.  The  ratio  of  peak  short-period  to  long-period  amplitude 
for  Bexar  is  0.379,  0.794  and  0.794  for  vertical,  radial  and  tangential  component,  respectively. 
For  Lee  Vining  the  ratio  of  peak  short-period  to  long-period  amplitude  is  0.221,  0.189  and 
0.0694  for  the  vertical,  radial  and  tangential  components.  For  Little  Skull  Mountain  the 
respective  ratios  are  0.0989,  0.394  and  0.455.  Again,  the  explosion  has  the  largest  short- 
period  to  long-period  amplitude  ratio  for  each  component. 

Another  comparative  difference  between  the  explosions  and  earthquakes  in  these  figures 
is  the  spectral  content  of  the  Rayleigh  waves.  The  explosion  Rayleigh  waves  exhibit  large, 
relatively  short-period  (3  to  8  sec.)  ’’ringing”  or  coda  waves  after  the  Airy  phase,  whereas 
the  earthquakes  display  primarily  the  dominant  long-period  Airy  phase.  This  effect  is  be- 
lievrd  to  be  a  dept  h  dependent  phenomena.  Kafka  (1990)  found,  in  a  study  of  New  England 
earthquakes  and  quarry  blasts,  that  shallow  events  produce  larger  short-period  (0.4  to  2.5 
II/.  /(*,,  waves  than  do  deeper  events.  Such  Rayleigh  wave  information  is  useful  as  a  depth 
dis<  riminant  to  he  p  distinguish  source  types. 

!  hose  M*t s  of  records  corroborate  I  lie  observat  ion  that  explosions  are  richer  in  short-period 
energy  relative  to  long  period  energy  as  compared  to  earthquakes.  Events  recorded  at  both 
FAS  and  CSC  ha  e  been  plotted  in  order  to  illustrate  that  this  observation  is  path  inde- 
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pendent.  Of  the  network  shown  in  Figure  1,  only  PAS  and  GSC  recorded  the  Lee  Vining 
event.  We  shall  show,  using  a  large  data  set,  that  this  spectral  difference  becomes  even  more 
evident  after  applying  propagational  corrections  and  can  be  used  as  an  effective  regional 
discriminant. 

Data  Analysis  and  Results 

We  compiled  seismic  moments  (;V/0)  and  local  magnitudes  (Ml)  for  NTS  explosions  and 
earthquakes  throughout  the  western  United  States  and  northern  Baja  California,  Mexico 
from  many  sources.  Mi  values  were  taken  from  the  CIT  and  Berkeley  catalogs,  and  from 
the  Northwestern  Mexico  Seismic  Network  (Vidal  and  Munguia,  1991).  For  recent  events 
(1988-1992)  in  Nevada  and  the  California-Nevada  Border  region,  Ml's  were  determined 
from  the  array  of  eight  broadband  stations  shown  in  Figure  1.  To  calculate  these  Aft’s,  an 
at  tenuation  curve  developed  by  Kanamori  el  al  (1992)  was  used.  All  the  Ml's  used  were 
determined  in  essentially  the  same  fashion. 

The  seismic  moments  collected  for  this  study,  however,  were  determined  through  a  variety 
of  means.  Dreger  and  Hefmberger  (1990.  1991a,  1991b,  &  1992),  Ma  and  Kanamori  (1991), 
and  Zhao  and  Heltnberger  (19911)  inverted  local  and  regional  broadband  waveforms  to  ob¬ 
tain  source  parameters,  including  moment.  Colin  ct  nl  (1982)  forward  modeled  near-field 
broadband  records  to  determine  earthquake  moments.  Moments  determined  fiom  combined 
regional  phase  (/’„/)  and  teleseismir  body  waves  (P  and  SB)  studies  were  also  gathered.  Bent 
and  Helmberger  (1991)  used  a  combination  of  forward  modeling  comparative  amplitude  ra¬ 
tios  l.o  estimate  moments  for  historical  southern  California  We  also  included  moments  and 
local  magnitudes  compiled  from  various  source  studies  for  historical  western  United  States 
earthquakes  (l)oser  and  Smith,  1989;  Doser,  1990). 
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Moments  from  '-everal  short-period  array  source  studies  were  incorporated  into  our  data 
set.  Mori  and  Firnkel  (1990)  obtained  moment  estimates  from  deconvolved  displacement 
pubi  s.  Employing  the  method  of  Brnne  (1970),  several  studies  calculated  Mo  from  the 
low  Irequency  amplitude  ol  shear  wave  spectra  (Johnson  and  McKvilly,  1974;  Fh'tcher  cl.  al, 
1984;  Frankel,  1984;  Vidal  and  Munguia,  1991). 

Surface  wave  moments  were  also  obtained  from  a  variety  of  sources.  Thio  and  Kanamori 
(1992)  obtained  source  parameters  from  broadband  TERRAscope  data  for  earthquakes 
throughout  southern  California  for  a  wide  range  of  magnitudes.  Their  inversion  method 
uses  both  Rayleigh  wave  and  Love  wave  spectra  and  employs  the  technique  of  Kanamori  and 
("liven  ( 1981).  Patton  and  Zandt  ( 1991 )  determined  moment  tensor  solutions  for  earthquakes 
t  hroughout.  the  western  U.  S.  using  a  linear  moment  inversion  scheme  developed  for  Rayleigh 
wave  spectra  by  Romanowicz  (1982)  and  extended  to  Love  wave  spectra  by  Patton  (1988). 
Wyss  and  Brune  (1968)  determined  moments  from  Love  wave  spectra  for  central  California 
events. 

Table  1  is  a  list  of  the  local  magnitudes  and  log-moments  for  earthquakes  for  which  source 
parameters  were  determined  from  broadband  data.  It  includes  all  available  data  from  the 
studies  by  Dreger  and  Helmberger  (1990,  1991a,  1991b,  1992),  Ma  and  Kanamori  (1991), 
and  I  hio  and  Kanamori  (1992). 

Table  2  lists  of  the  local  magnitudes  and  log-moments  for  explosions  cited  in  this  study. 
Seven  explosion  moments  were  determined  from  near-field  observations  (Aki  et  al,  1974; 
Helmberger  and  Hadley.  1981;  Stump  and  Johnson,  1984;  and  Johnson,  1988),  although 
most  were  obtained  from  surface  wave  studies.  Stevens  (1986)  calculated  spectral  moments 
from  Rayleigh  wax e  spectra  for  large  NTS  explosions.  Given  and  Mellman  (1986)  performed 
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moment  tensor  inversions  of  large  NTS  blasts,  too,  using  the  path  structures  developed  by 
Stevens  (1986).  They  used  Rayleigh  wave  and  Love  wave  spectra  to  solve  for  the  isotropic 
source  as  well  as  a  double  couple  source  associated  with  tectonic  release.  The  data  sets  for 
these  last  two  studies  overlap  substantially.  The  moment  values  were  found  to  be  similar, 
so  that  for  events  for  which  two  moments  were  available,  the  Stevens’  ( 1 986 J  moment  was 
used. 

Time  domain  surface  wave  moments  of  explosions  were  also  included  in  this  study.  Woods 
and  Harkrider  (1993)  determined  moments  from  the  peak  to  peak  amplitude  (PPA)  of  the 
dominant  Airy  phase  of  the  Rayleigh  wave  for  NTS  explosions.  With  their  technique,  moment 
is  determined  from  the  ratio  of  the  observed  PPA  to  that  of  a  synthetic  seismogram  with  a 
given  input  moment.  Figure  3  plots  comparisons  of  observed  vs.  modeled  vertical  Rayleigh 
waves  for  the  55  station  network  used  in  their  study.  Data  from  representative  Pahute 
Mesa  NTS  explosions,  with  little  or  no  tectonic  release,  are  the  upper,  darker  traces.  The 
synthetics  seismograms  were  generated  with  a  step  moment  source  buried  at  0.6  km,  a  typical 
shot  depth.  Many  of  the  path  structures  used  were  taken  from  the  Stevens  (1986)  study. 
Other  paths  were  determined  by  inverting  Rayleigh  wave  dispersion  data.  It  can  be  seen  that 
the  waveform  fits  are  good,  in  that  dispersion  and  amplitude  are  both  well-modeled.  The 
correlation  between  moments  obtained  by  this  time  domain  moment  and  a  more  standard 
spectral  scheme  performed  on  a  108  event  subset  of  this  study  is  very  good  (s<*  Woods  and 
Harkrider,  1993),  implying  that  this  time  domain  moment  method  yields  accurate,  robust 
moment  measurements.  Using  this  method  we  obtained  moments  for  other  small  NTS  events. 

Figure  4a  displays  Mr,  vs.  log  Moment  (N-m)  plotted  for  299  earthquakes  and  178  explo¬ 
sions.  Solid  symbols  represent  earthquakes,  while  open  symbols  and  crosses  represent  NTS 
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explosions.  Moments  determined  from  near-field  body  wave  studies  are  plotted  as  circles. 
Earthquake  moments  determined  from  surface  wave  studies  and  body  wave  studies  are  repre¬ 
sented  by  triangles  and  diamonds,  respectively.  Explosion  moments  determined  from  surface 
source  studies  are  denoted  by  stars  (Given  and  Mellman,  1986),  squares  (Stevens,  1986)  and 
crosses  (Woods  and  Harkrider,  1993).  This  figure  shows  how  well  this  discriminant  works. 
There  is  a  significant  separation  of  earthquakes  and  explosions,  with  no  real  overlap  of  the 
two  populations.  This  discriminant  also  works  at  all  scales,  with  explosions  and  earthquakes 
following  their  respective  scaling  laws  over  a  wide  range  of  magnitudes  and  moments;  for 
earthquakes  this  is  true  over  seven  and  a  half  orders  of  magnitude.  It  should  be  noted  that 
i  he  earthquakes  with  a  log(.V/„)  below  13.0  were  determined  from  local  stations  (D  <  75  km) 
and  would  not  be  detectable  at  regional  distances.  They  are  included  here  only  to  show  the 
continuity  of  the  linear  scaling  relationship  between  Ml  and  log(M0)  for  earthquakes. 

Figure  4b  is  a  blow-up  of  the  portion  of  Figure  4a  containing  explosion  data.  It  is  im¬ 
portant  to  note  when  examining  these  two  plots  that  the  data  are  taken  from  a  number  of 
sources.  M l ’s  were  determined  from  different  networks  and  the  moments  were  calculated  in 
a  valid v  of  ways.  Vet  a  distinct  separation  of  the  two  populations  is  still  obtained.  There  is 
one  anomalous  explosion.  Buggy  ( Ml  —  3.96,  log(Mo)=  14.28),  which  lies  very  close  to  the 
earthquake  population.  This  was  a  Plow  Shares  event  in  which  four  nuclear  charges  with  an 
announced  combined  yield  of  5.4  kt  were  detonated.  We  would  expect  that  a  multiple-source 
event  like  this  one  would  display  characteristics  similar  to  a  distributed  earthquake  source, 
he.  it  should  be  enriched  in  long-period  energy  relative  to  a  point  source  explosion  of  the 
same  yield. 
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Discussion  and  Conclusions 


The  Ml: Mo  criterion  appears  to  be  a  robust  method  to  discriminant  regional  seismic 
events.  For  a  given  moment,  the  Ml  of  an  explosion  is  more  than  0.5  units  larger  than 
that  of  an  earthquake.  This  difference  can  be  attributed  to  Ml  being  a  short-period  (1  Hz) 
energy  measurement,  whereas  the  moment  is  determined  from  long-period  body  wave  phases 
(period  >  4  seconds)  and  Rayleigh  waves  (10  to  40  seconds). 

There  are  several  explanations  for  the  observed  difference  in  short-period:long-period  spec¬ 
tral  character  between  earthquakes  and  explosions.  Earthquakes  tend  to  be  asperity-driven 
distributed  sources,  as  opposed  to  explosions  which  can  more  aptly  be  modeled  as  point- 
sources  with  impulsive  time  functions.  Dreger  and  Helmberger  (1991a)  showed  that  broad¬ 
band  seismograms  from  small  local  earthquakes  (4.0  <  Ml  <  5.0)  cam  be  modeled  as  dis¬ 
tributed  finite  sources  rather  than  as  point- sources.  Such  distributed  slip  time  functions  will 
generally  result  in  reduced  high-frequency  spectra  relative  to  a  point-source  step  moment. 
It  is  possible  that  even  very  small  earthquakes  (Ml  <4)  behave  similarly.  Another  possible 
reason  that  this  discriminant  works  at  low  magnitudes  is  that  for  small  events  the  quantities 
being  measured  are  the  P-wave  low  frequency  spectrum  (Ml)  and  the  S-wave  low  frequency 
spectrum  (Mo),  in  which  case,  for  explosions,  the  S-wave  energy  is  predominantly  generated 
by  secondary  P  to  S  converted  phases.  If  this  is  so,  the  moment  should  not  be  determined 
solely  from  P-wave  information. 

Spall  effects  may  also  be  a  cause  of  observed  high-frequency  (0.15  Hz  <  /  <  2.0  Hz) 
enrichment  of  explosion  source  spectra.  Theoretical  results  (Day  and  McLaughlin,  1991) 
and  various  observational  studies  ( Viccclli,  1973;  Stump,  1985;  Taylor  and  Randall,  1989; 
and  Patton,  1990)  conclude  that  spallation  can  be  a  significant  contributor  to  short  period 
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energy,  while  Patlon  (1988)  and  Day  ct  al  (1983)  find  no  appreciable  spall  energy  at  periods 
greater  than  8  sec  onds.  S pall  energy  would  tend  to  increase  Ml  measurements. 

'This  discriminant  is  only  limited  bv  t  he  detection  threshold  capability  of  long-period  data, 
as  t  he1  two  populat  ions  do  not  converge  at  small  magnitudes.  This  observation  implies  that 
the  convergence  in  the*  Ms'-m-b  ratio  for  small  magnitude  events  seen  in  some  previous  studies 
is  due  to  approaching  the  effective  signal  to  noise  level  for  measuring  surface  amplitudes. 
Previous  studies  relied  on  data  recorded  on  lower  grade  (usually  analog)  instruments  that 
did  not  have  the  recording  capabilities  of  modern,  high  dynamic  range,  digital  seismometer 
systems.  Modern  data  combined  with  digital  processing  techniques  increases  the  resolution 
of  long-period  transient  signals.  Figure  5  shows  the  Yucca  Flat  blast  Floydada  recorded  at 
four  TERRAscope  stations  (epicentral  distances  being  between  200  and  400  km),  played-out 
with  PE  and  WA  instruments.  The  amplitudes  indicate  that  these  signals  would  not  be 
discernible  on  the  actual  analog  instruments.  Ml  for  this  event  is  4.0  and  its  log  moment  is 
14.20.  Assuming  it  is  a  shallow  explosion  above  the  water  table,  the  yield  can  be  inferred 
to  be  less  than  10  kt  from  the  moment-yield  scaling  relationships  determined  for  NTS  by 
Woods  and  Harkrider  (1993).  Were  it  detonated  in  hard  rock  below  the  water  table,  it  would 
correspond  to  a  two  kiloton  explosion.  We  estimate  that  were  this  event  2.5  times  lower  in 
yield  it  would  still  be  possible  to  obtain  its  moment,  yielding  a  magnitude  threshold  of  Ml  > 
3.6  for  explosions. 

Low  SNR  seismograms  also  may  be  phase  match  filtered  (Herrin  and  Goforth,  1977;  and 
Stevens,  1986a)  in  order  to  retrieve  the  signal  of  very  small  events  for  spectral  amplitude 
estimates.  However,  there  is  some  debate  on  the  accuracy  of  using  this  method,  see  Der 
(1986)  and  Stevens  (1986b)  for  a  discussion  of  this  problem.  Employing  one  or  another 
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means  of  time  domain  measurements  in  order  to  obtain  the  seismic  moment  avoids  such 
problems. 

Simple  time  domain  moment  measurements  included  in  this  study  are  straightforward 
and  can  be  directly  applied  to  historical  analog  data  sets  in  order  to  establish  earthquake 
Ml'.Mq  curves  for  other  regions.  Since  only  peak  amplitude  measurements  are  necessary 
for  such  methods,  the  required  signal  to  noise  ratio  (SNR)  effectively  decreases,  too.  Time 
domain  peak  amplitude  measurements  are  also  less  susceptible  to  noise  contamination  than 
are  spectral  amplitude  measurements.  The  other  parameter  Ml  is  a  simple  time-domain 
measurement  which  can  be  made  on  the  smallest  detectable  seismic  sources.  Modem  source 
inversion  techniques  that  make  use  of  regional  body  wave  phases  (such  as  P^  and  Snj) 
recorded  at  sparse  regional  networks  also  make  for  powerful  moment  determination  tools  as 
shown  in  studies  using  the  TERRAscope  broadband  network  (Dreger  and  Helmberger,  1992; 
and  Zhao  and  Helmberger,  1993). 

A  disadvantage  to  determining  moment  tensor  solutions  is  that  many  methods  require 
(Jreen’s  functions  in  order  to  obtain  accurate  moments.  However,  with  modern  broadband, 
high  dynamic  range  instrumentation  it  is  quite  feasible  to  use  moderate-sized  to  small  events 
to  determine  regional  path  structures.  Zhao  and  Helmberger  (1991)  detail  the  forward  model¬ 
ing  of  Pni ,  Sni  and  Rayleigh  wave  regional  phases  along  a  continental  shield  path.  Dreger  and 
Helmberger  (1990)  were  able  to  forward  model  velocity  structure  using  waveform  data  from 
small  (mj,  =  3.7  -  4.0)  local  events.  Shallow  crustal  structure  can  also  be  inferred  from  the 
inversion  of  surface  wave  dispersion  data  generated  by  small,  regional  seismic  sources  (Saikia 
cl  at.  1990)  as  well  as  in  conjunction  with  teleseismic  surface  waves  (Thio  and  Kanamori, 
1991).  These  path  structures  needn’t  be  overly  complicated  or  detailed.  Dreger  and  Helm- 
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berger  (1991b,  1992)  found  that  they  could  model  broadband  regional  body- wave  phases 
from  events  located  throughout  Southern  California  and  recorded  at  various  TERRAscope 
stations  with  one  relatively  simple  crustal  model  (with  an  underlying  mantle  half-space). 
Using  their  regional  (I reerrs  functions  iri  conjunction  with  a  time-domain  source  inversion 
method  yielded  source  parameter  solutions  consistent  with  other  studies.  Zhao  and  Helm- 
berger  (1993)  extended  this  technique  to  include  surface  waves  and  found  that  this  inversion 
scheme  worked  well  with  simple  regional  eart  h  models  as  well. 

An  important  point  to  address  concerns  the  a  priori  source  type  assumptions  made  in 
determining  the  seismic  moments.  All  earthquake  moments  were  determined  assuming  ei¬ 
ther  a  double  couple  source  (for  moment  tensor  and  Green’s  function  inversion  methods) 
or  a  circular  fault  model  (for  corner  frequency  moment  estimates).  Most  of  the  the  NTS 
explosion  moments  were'  determined  by  modeling  the  Rayleigh  waves  as  being  generated  by 
a  shallow,  isotropic  source  excitation  function  in  the  absence  of  an  azimuthally  dependent 
radiation  pattern.  The  Given  and  Mellman  (1986)  moments,  however,  were  determined  by 
inverting  Rayleigh  and  Love  wave  spectra  to  obtain  the  isotropic  and  double  couple  com¬ 
ponents  of  the  moment  tensor  solution.  An  obvious  problem  with  classifying  events  with 
respect  to  their  Mf.Mo  ratio  is  that  for  an  unidentified  seismic  event  one  doesn’t  know 
which  source  model  to  assume  in  order  to  estimate  the  moment.  This  point  leads  to  the 
question,  how  distinguishable  would  the  two  populations  be,  had  they  all  been  treated  as 
double  couple  sources  for  moment  estimation  purposes?  For  a  given  seismic  moment,  the  av- 
erage  radiation  pat  tern  amplitude  (assuming  absolute  values)  for  Rayleigh  waves  generated 
by  a  predominantly  strike-slip  fault  motion  (sin20  radiation  pattern)  is  only  slightly  larger 
(by  10  percent)  than  that  of  an  explosion.  In  this  case  one  could  still  infer  the  source  type 
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from  the  moment  estimate.  For  such  strike-slip  earthquakes,  observations  from  all  azimuthal 
quadrants  (the  four  lobes  of  the  .>in20  radiation  pattern)  will  result  in  better  constrained 
moment  determinations. 

'The  difficulty  occurs  in  modeling  the  Rayleigh  waves  assuming  a  shallow  source  with  a 
dip-slip  orientation  for  which  there  is  a  singularity  in  the  Green’s  function  solutions.  In 
this  special  case  we  would  obtain  a  larger  moment  and  the  event’s  Mi  -  M0  ratio  would 
decrease  towards  the  earthquake  population.  However,  one  can  generally  model  and  invert 
for  earthquake  sources  that  occur  deeper  than  3  km.  When  an  event  does  not  fit  our  regional 
Green’s  functions  to  some  specified  degree,  we  would  assume  it  is  shallow  and  model  it  as  a 
vertical  strike-slip  earthquake.  Thus,  explosions  would  still  discriminate  as  displayed  above 
and  shallow  strike-slip  earthquakes  would  still  lie  within  in  the  earthquake  population.  M0 
for  a  shallow  dip-slip  earthquake,  however,  may  be  underestimated  and  could  potentially  fall 
within  the  explosion  population.  In  a  monitoring  environment  we  would  have  a  problematic 
event,  to  which  other  discrimination  criterion  would  need  to  be  applied  in  order  to  identify 
it  correctly. 

Applying  the  Mi'.Mq  ratio  in  conjunction  with  other  discriminants  would  yield  a  more 
effective  source  identification  scheme.  Phase  information  (i.e.  Rayleigh  wave  polarity)  and 
Love  wave  data  would  be  helpful  in  such  cases,  since  reversed  polarity  Rayleigh  waves  and/or 
large  Love  wave  amplitudes  are  diagnostic  of  a  double  couple  source.  A  depth  discriminant, 
based  on  the  complexity  and  amplitude  of  the  Rayleigh  wave  tail  (coda)  for  example,  would 
be  useful  for  depth  constraints  lo  be  placed  on  the  moment  tensor  inversion  of  a  source. 
In  a  related  study  Woods  and  Ilelmberger  (1992)  found  that  the  ratio  of  sln>rt-»  eriod  en¬ 
ergy  in  the  vertical  component  F’nl  wave  train  to  that  in  the  long-period  surface  wave  train 
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(siimmcd  over  all  three  components)  of  regional  seismograms  also  separates  source  popula¬ 
tions,  although  it,  too,  suffers  similar  magnitude  threshold  limitations. 

For  very  small  events  information  from  high-frequency  (f  >  1  Hz)  phases  may  be  needed 
to  compliment  the  A//,:A/0  discriminant.  Comparisons  of  Lg  spectral  amplitude  levels  for 
different  bandwidths  (in  the  0.5  to  8  11/  range)  have  been  shown  to  effectively  discriminate 
events  with  magnitudes  down  to  m/,,/.n)  =  A. A  (Murphy  and  Bennett,  1982);  Taylor  et  al, 
1 9SS >.  Although  the  Mi,:Mu  discrimination  threshold  is  slightly  higher  than  this,  the  regional 
body  wave  phases  and  Rayleigh  waves  used  in  this  method  are  not  as  susceptible  to  path 
"blockage'’  effects  as  is  the  phase.  Also,  using  relatively  close-in  stations  for  surface 
wave  measurements,  as  it  is  possible  to  do  with  this  method,  path  attenuation  effects  are 
minimized.  The  Mi/.Mo  ratio  discriminant  would  work  well  as  a  companion  test,  or  check, 
for  other,  high-frequency,  discrimination  methods,  particularly  for  events  in  the  3.5  <  m*  < 
4.5  range  for  which  teleseismic  methods  no  longer  work. 

Because  both  source  parameters  used  in  this  discrimination  method  can  be  obtained  from 
a  sparse  broadband  network,  this  discriminant  can  be  applied  throughout  the  world  as  more 
broadband  stations  similar  to  those  of  the  IRIS  network  come  on  line.  For  an  active  tectonic 
region  the  threshold  for  this  discriminant  is  A//.  =  A.l  for  earthquakes  and  Mi  —  3.6  for 
explosions  for  epicentral  distances  up  to  600  km. 
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Figure  1.  Map  showing  the  study  area.  Broadband  stations  used  to  determine  Mq  and 
Ml  for  1988+  events  are  shown  with  triangles.  Events  shown  in  this  study  are  denot¬ 
ed  by  stars. 

Figure  2.  (a)  Broadband  displacement  (middle  trace)  records  of  the  JVE  event  Kear- 
sarge  and  the  Lee  Vining  and  Skull  Mountain  earthquakes  recorded  at  PAS  and  played 
out  with  long-period  (top  trace)  and  short-period  Gower  trace)  instruments,  (b)  Analo¬ 
gous  plot  for  the  explosion  Bexar  and  the  Lee  Vining  and  Skull  Mountain  earthquakes 
recorded  at  GSC. 

Figure  3.  Comparison  of  synthetic  and  observed  fundamental  Rayleigh  waves  for  a 
regional  network  from  which  seismic  moments  were  determined.  The  data  arc  the 
upper,  darker  traces.  The  records  chosen  came  from  an  ensemble  of  high-SNR  Pahute 
Mesa  events  with  little  or  no  tectonic  release.  The  synthetics  were  generated  with  a 
step  moment  isotropic  source  buried  at  0.6  km. 

Figure  4.  (a)  ML  vs.  log  M0  for  299  earthquakes  (solid  symbols)  and  178  explosions 
(open  symbols  and  crosses).  Circles  arc  moments  from  near-field  source  studies.  Tri¬ 
angles  are  surface  wave  moments  and  diamonds  are  body  wave  moments.  Explosion 
moments  determined  from  surface  waves  are  shown  as  stars  (Given  and  Mellman, 
1985),  squares  (Stevens,  1986a)  and  crosses  (Woods  and  Harkrider,  1992).  (b)  Blow¬ 
up  of  that  part  of  the  plot  which  contains  explosions. 
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Figure  5.  A  small  Yucca  Flat  explosion,  Floydada  (ML=4.0),  recorded  at  TERRA- 
scope  stations  and  played  out  cm  Press-Ewing  (long-period)  and  Wood-Anderson 
(short-period)  instmments. 
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Floydada  (Yucca  Flats),  Ml=4.0 
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Table  1 :  Earthquakes  with  Moments  Determined  Using  Broadband  Data 


Date 

log(M0) 

Event 

Date 

Ml 

log(M0) 

Event 

88/06/10 

5.4 

17.11  ° 

Gorman 

92/03/03 

14.49  ° 

Walker  Pass 

88/06/26 

4.6 

15.78  d 

Chino 

92/03/04 

4.2 

15.14  ° 

San  Clemente 

88/06/27 

5.5 

16.77  c 

San  Juan  Botista 

92/03/05 

3.8 

14.59° 

Bakersfield 

88/07/06 

3.7 

14.88  d 

Chino  A.S. 

92/04/10 

3.4 

14.60° 

Borrego 

88/12/03 

4.9 

16.38  h 

Pasadena  M.S. 

92/04/15 

3.4 

14.00“ 

Lytle  Creek 

11:40: 

2.1 

12.41  h 

Pasadena  A.S.  1 

92/04/23 

4.6 

15.49  ° 

Joshua  Tree  F.S. 

12:08: 

2.0 

13.02  h 

Pasadena  A.S.  3 

92/04/23 

6.1 

18.29° 

Joshua  Tree  M.S. 

12:13: 

2.4 

13.02  h 

Pasadena  A.S.  4 

13:58: 

4.1 

15.01  ° 

II 

It 

A.S.  1 

12:15: 

1.6 

12.40  b 

Pasadena  A.S.  5 

22:55: 

3.8 

14.54  ° 

M 

II 

A.S.  2 

13:36: 

1.8 

12.09  h 

Pasadena  A.S.  6 

23:52: 

3.8 

14.54“ 

It 

II 

A.S.  3 

14:46: 

1.9 

11.88  h 

Pasadena  A.S.  7 

92/04/24 

3.5 

14.37  ° 

II 

II 

A.S.  4 

88/12/04 

2.0 

12.24  h 

Pasadena  A.S.  8 

18:06: 

3.7 

14.52  ° 

M 

II 

A.S.  5 

88/12/08 

2.2 

12.64  b 

Pasadena  A.S.  9 

18:20: 

3.7 

14.61  ° 

It 

II 

A.S.  6 

89/01/19 

5.0 

16.51  ° 

Malibu 

92/04/25 

3.7 

14.73  ° 

H 

II 

A.S.  7 

89/02/18 

4.3 

15.34  f 

Upland 

18:56: 

4.4 

15.05  ° 

l» 

II 

A.S.  8 

89/08/08 

5.3 

16.41  c 

Los  Gatos 

92/04/26 

3.7 

14.27  ° 

II 

II 

A.S.  9 

89/10/18 

7.0 

19.48  c 

Loma  Prieta 

6:26: 

4.2 

15.69° 

II 

II 

A.S.  10 

90/02/28 

3.7 

14.65  ' 

Upland  F.S. 

17:21: 

4.3 

15.39° 

H 

II 

A.S.  11 

90/02/28 

5.2 

17.40  ' 

Upland  M.S. 

92/04/27 

4.2 

15.41  ° 

II 

It 

A.S.  12 

90/03/01 

4.7 

15.70  f 

Upland  A.S. 

92/04/28 

3.7 

14.90° 

II 

II 

A.S.  13 

90/03/02 

4.6 

15.60  ' 

Upland  A.S. 

11:33: 

3.8 

15.12  ° 

II 

II 

A.S.  14 

90/04/17 

4.6 

15.78  f 

Upland  A.S. 

92/04/30 

B 

14.70° 

II 

II 

A.S.  15 

Table  1  (corn.):  Earthquakes  with  Moments  Determined  Using  Broadband  Data 

Di 

Ml 

log  (M0) 

Event 

Date 

Ml 

log(M0) 

Event 

9(V10/24 

5.0 

16.70  * 

Lee  Vining 

92/05/01 

3.8 

14.71  « 

J.  T.  A.S.  16 

90/12/17 

3.7 

14.38  ° 

Big  Bear  I 

92/05/02 

4.1 

14.64° 

"  "  A.S.  17 

90/12/18 

3.7 

15.04  * 

White  Wolf  Fault 

19:10: 

3.4 

14.35  ° 

"  "  A.S.  18 

91/05/20 

3.7 

14.48  ° 

San  Jacinto  I 

92/05/04 

4.0 

15.04° 

"  "  A.S.  19 

91/05/20 

3.7 

14.20  ° 

San  Jacinto  II 

16:19: 

4.8 

16.19° 

"  ”  A.S.  20 

91/06/28 

5.4 

17.41  f 

Sierra  Madre 

92/05/06 

15.87  “ 

"  "  A.S.  21 

15:37: 

3.9 

14.61  / 

S.  Madre  A.S. 

92/05/12 

m 

15.60° 

"  "  A.S.  22 

4.3 

15.60  ? 

S.  Madre  A.S. 

92/05/18 

3.5 

14.39  ° 

"  "  A.S.  23 

91/06/29 

4.0 

14.60° 

Mojave 

15:44: 

4.9 

16.10° 

"  "  A.S.  24 

91/07/06 

3.8 

14.59  f 

S.  Madre  A.S. 

92/06/11 

4.3 

15.29  ° 

"  "  A.S.  25 

91/10/12 

4.0 

15.04° 

Blue  Cut 

92/05/31 

3.2 

13.78  ° 

Lenwood  Flat  1 

91/10/27 

14.26° 

San  Jacinto 

92/05/31 

3.5 

14.28  ° 

Lcnwocxl  Flat  11 

91/12/03 

16.72  * 

San  Miguel 

92/06/28 

3.7 

14.29  ° 

Landers  A.S. 

91/12/04 

14.95  ° 

Julian 

92/06/29 

5.4 

16.85  ° 

Landers  A.S. 

91/12/04 

4.0 

14.62  ° 

Big  Bear  II 

14:41: 

4.4 

16.90° 

Landers  A.S. 

92/02/17 

3.5 

14.66  ° 

Coso 

16:01: 

5.2 

17.16° 

Landers  A.S. 

92/02/19 

4.0 

14.85  ° 

Coso  I 

92/06/29 

5.6 

17.40  c 

Little  Skull  Mtn 

92/02/19 

3.7 

14.41  ° 

Coso  II 

92/06/29 

3.7 

14.66“ 

La  Canada 

92/02/21 

3.7 

15.30  ° 

Coso 

92/06/30 

4.7 

15.83  * 

Landers  A.S. 

92/02/22 

3.9 

14.28  " 

Coso 

21:49: 

4.3 

15.87  ° 

Landers  A.S. 

Thio  and  Kanamori,  1992;  Ma  and  Kana:  ori,  1991; (c)  Zhao  and  Helmberger,  1993: 

(</).  ro.  t/).  <g )  Dreger  and  Helmberger,  1990,  1991a,  1991b,  1992; 

36 


Table  2:  Explosions  in  this  study 


Date 

*4 

log(M0) 

Event 

Date 

A4 

log(Af0) 

Event 

57/09/19 

4.25 

14.30° 

Rainier 

70/03/19 

4.05 

13.65  * 

Jal 

61/12/03 

3.91 

13.86° 

Fisher 

70/03/23 

5.29 

15.73  * 

Shaper 

62/06/27 

4.63 

14.55  * 

Haymaker 

70/03/26 

6.17 

17.72  c 

Handley 

65/02/16 

4.21 

13.98  * 

Merlin 

70/05/01 

4.18 

14.35  h 

Beeblam 

65/03/03 

5.04 

15.41  * 

Wagtail 

70/05/05 

4.82 

14.70  * 

Mintleaf 

65/04/14 

4.23 

14.02  * 

Palanquin 

70/05/15 

5.00 

15.19* 

Cornice 

65/05/12 

3.91 

13.84  * 

Buteo 

70/05/21 

3.60 

13.23  * 

Manzanas 

66/04/25 

4.69 

14.90  * 

Pinstripe 

70/05/21 

5.00 

14.98  * 

Morrones 

66/04' 14 

5.02 

15.39  * 

Duryea 

70/05/26 

4.63 

14.60* 

Hudsonmoon 

66/05/05 

4.04 

14.24* 

Cyclamen 

70/05/26 

5.20 

15.56  * 

Flask 

66/05/13 

5.01 

15.75  * 

Piranha 

71/06/19 

4.20 

13.83  * 

Embudo 

66/06/02 

4.98 

15.73  * 

Piledriver 

71/06/23 

4.50 

14.45  * 

Laguna 

66/06/25 

4.33 

14.13  * 

Vulcan 

71/06/29 

4.90 

14.64  * 

Harebell 

66/06/30 

5.90 

16.49  * 

Halfbeak 

71/10/29 

4.10 

14.03  * 

Pedemal 

66/12/20 

6.07 

16.98  * 

Greeley 

71/10/08 

4.10 

14.02  * 

Cathay 

67/02/08 

4.09 

13.82  * 

Ward 

72/04/19 

4.20 

13.73  * 

Longchamps 

67/02/23 

3.89 

13.81  * 

Persimmon 

72/05/19 

4.48 

14.66* 

Monero 

67/04/07 

3.97 

13.76  * 

Fawn 

72/10/03 

4.40 

13.93  * 

Delphinium 

67/04/21 

3.83 

13.51  h 

Effendi 

73/06/06 

5.68 

16.78  ' 

Almendro 

67/04/21 

3.89 

13.72  h 

Chocolate 

75/05/14 

5.86 

16.50  « 

Tybo 

67/05/23 

5.57 

16.46  * 

Scotch 

75/06/19 

5.74 

16.58  * 

Mast 

67/06'26 

4.54 

14.42  * 

Midimist 

75/06/26 

5.95 

16.96  * 

Camembert 
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Table  2  (cont.):  Explosions  in  this  study 


Date 

A4 

log(M0) 

Event 

Date 

*4 

log(Af  0) 

Event 

67/08/31 

4.73 

14.46  * 

Doormist 

75/10/24 

4.86 

14.80* 

Huskypup 

67/09/21 

3.87 

13.51  * 

Marvel 

75/10/28 

6.08 

17.03  * 

Kasseri 

67/09/27 

5.66 

16.32  ' 

Zaza 

76/02/12 

6.12 

17.00* 

Fontina 

68/01/19 

5.74 

16.83  * 

Faultless 

76/02/26 

5.86 

16.79  e 

Cheshire 

68/03/12 

3.96 

14.28  * 

Buggy 

76/03/09 

5.83 

16.79  ' 

Estuary 

68/04/26 

6.09 

17.24  * 

Boxcar 

76/03/17 

5.77 

16.52  ' 

Pool 

68/09/24 

5.06 

14.85  * 

Hudsonseal 

77/04/05 

5.20 

16.01  e 

Marsilly 

68/11/04 

4.46 

14.57  * 

Crew 

77/04/27 

4.90 

15.62  * 

Bulkhead 

68/12/19 

6.14 

17.35  * 

Benham 

77/05/25 

5.00 

15.52  * 

Crewline 

69/02/12 

4.81 

14.63  * 

Cypress 

77/11/09 

5.40 

16.27  * 

Sandreef 

69/05/07 

5.68 

16.42  ' 

Purse 

77/12/14 

5.30 

16.22  e 

Farallones 

69/09/16 

6.19 

16.82  b 

Jorum 

78/03/23 

5.30 

16.02  d 

Iceberg 

69/10/08 

5.50 

16.19  c 

Pipkin 

78/04/11 

5.10 

15.89  d 

Backbeach 

69/10/29 

4.40 

13.55  * 

Cruet 

78/07/12 

5.20 

15.75  * 

Lowball 

69/10/29 

4.60 

13.97  * 

Pod 

78/08/31 

5.30 

15.92  d 

Panir 

69/10/29 

5.50 

15.79  * 

Calabash 

78/09/13 

4.40 

14.55  * 

Diablohawk 

70/02/11 

4.67 

14.77  * 

Diamamist 

78/09/27 

5.30 

16.15  d 

P.ummy 

70/02/25 

4.60 

14.89  * 

Cumarin 

78/11/02 

4.23 

14.36  h 

Emmenthal 

70/02/26 

4.80 

14.65  * 

Yannigan 

78/11/18 

5.00 

15.30  * 

Quargel 

70/03/06 

4.20 

13.57  * 

Cyathus 

78/12/16 

5.20 

16.40  c 

Farm 

70/03/06 

4.10 

13.09  * 

Arabis 

79/02/08 

5.20 

15.80* 

Quinella 

79/06/1 1 

5.30 

16.13  d 

Pepato 

84/03/31 

4.20 

13.71  * 

Agrini 

Table  2  (cont.):  Explosions  in  this  study 


Date 

A4 

log(M0) 

Event 

Date 

Ml 

log(M0) 

Event 

79/06/20 

3.92 

13.91  h 

84/05/01 

5.10 

15.76  * 

Mundo 

79/09/06 

5.30 

16.16  d 

84/05/31 

5.30 

15.89  * 

Caprock 

79/09/26 

5.20 

15.98  d 

Sheepshead 

84/06/20 

4.40 

14.48  * 

Duoro 

79/11/29 

3.89 

13.66  h 

Backgammon 

84/07/25 

5.30 

15.42  * 

Kappeli 

80/03/08 

3.89 

13.43  * 

Norbo 

84/08/02 

4.30 

14.28  * 

Correo 

80/04/26 

5.10 

15.91  f 

Col  wick 

84/08/30 

4.50 

14.52  * 

Dolcetto 

80/04/30 

5.10 

15.61  * 

Pyramid 

84/09/13 

4.80 

15.01  * 

Breton 

80/07/25 

5.10 

16.10  d 

Tafi 

84/11/10 

4.30 

14.21  * 

Villita 

80/10/31 

4.50 

14.77  * 

Minersiron 

84/12/09 

5.10 

15.63  * 

Egmont 

81/01/15 

5.20 

15.79* 

Baseball 

84/12/15 

5.00 

15.63  * 

Tierra 

81/06/06 

5.40 

15.87/ 

Harzer 

85/03/15 

4.60 

14.77  * 

Vaughn 

81/11/12 

5.00 

15.61  * 

Rousanne 

85/03/23 

5.00 

15.31  * 

Cottage 

82/01/28 

5.30 

16.08  * 

Jornada 

85/04/02 

5.40 

16.12  * 

Hermosa 

82/02/12 

5.00 

15.95  * 

Molbo 

85/04/06 

4.50 

14.86  * 

Misty-rain 

82/02/12 

5.00 

15.90* 

Hosta 

85/06/12 

5.10 

15.85  * 

Salut 

82/04/17 

4.40 

14.35  * 

Tenaja 

85/07/25 

5.20 

15.73  * 

Serena 

82/04/25 

5.00 

15.79* 

Gibne 

85/12/05 

5.20 

15.56* 

Kinibito 

82/05/07 

5.20 

15.76* 

Bouschet 

85/12/28 

5.00 

15.57  * 

Goldstone 

82/06/24 

5.20 

15.96  * 

Nebbiolo 

86/03/22 

5.10 

15.46* 

Glencoe 

82/07/29 

4.10 

14.45  * 

Monterey 

86/04/10 

4.50 

14.78  * 

Mighty  oak 

82/08/05 

5.40 

16.16  * 

Atrisco 

86/04/22 

5.00 

15.71  * 

Jefferson 

82/09/02 

3.50 

13.23  * 

Cerro 

86/05/21 

4.00 

13.96* 

Panamint 
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Table  2  (corn.):  Explosions  in  this  study 


Date 

Ml 

log(Afo) 

Event 

Date 

Ml 

log(A/0> 

Event 

82/09/23 

4.60 

14.82  « 

Huronlanding 

86/06/05 

5.30 

15.58* 

Tajo 

82/09/23 

4.60 

14.88  * 

Frisco 

86/06/25 

5.30 

15.61  * 

Darwin 

82/09/29 

3.80 

13.98  * 

Boirego 

86/07/17 

5.20 

15.08* 

Cybar 

82/11/12 

4.10 

13.92  * 

Seyval 

86/07/24 

4.70 

14.40* 

Cornucopia 

82/12/10 

4.40 

14.56  * 

Manteca 

86/09/30 

5.40 

15.72  * 

Labquark 

83/03/26 

4.90 

15.44* 

Cabra 

86/10/16 

5.30 

15.77  * 

Belmont 

83/04/14 

5.20 

15.59  * 

Torquoise 

87/04/18 

5.30 

15.77  < 

Delamar 

83/05/26 

4.30 

14.26* 

Fahada 

87/04/30 

5.30 

15.88  * 

Hardin 

83/06/09 

4.30 

14.29  * 

Danablu 

87/06/18 

4.00 

13.90  h 

Brie 

83/09/01 

5.30 

15.77/ 

Chancellor 

87/10/23 

5.40 

15.31  * 

Borate 

83/09/22 

4.00 

13.90* 

Techado 

88/04/07 

3.60 

12.75  * 

Abliene 

83/12/16 

4.80 

15.22  * 

Romano 

90/06/21 

4.30 

13.91  h 

Austin 

83/05/05 

4.20 

14.26  * 

Crowdie 

91/08/15 

4.00 

14.10  h 

Floydada 

84/02/15 

4.50 

14.78  * 

Milagro 

92/09/18 

4.00 

14.02  h 

Hunterstrophy 

84/03/01 

5.30 

15.85  * 

Tortugas 

92/09/23 

4.41 

14.27  * 

Divider 

(a)  Aki  et  al,  1974;  Helmberger  and  Hadley,  1981;  (c)  Stump  and  Johnson,  1984; 

(</)  Given  and  Mellman,  1986;  (e)  Stevens,  1986a;  Johnson,  1988; 

<*>  Woods  and  Harkrider,  1993;  (fc)  This  study 
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Pasadena,  CA  91125,  USA 


ABSTRACT 

A  method  of  relocation  and  source  characterization  of  small  earthquakes 
using  one  modern  regional  station  is  developed  and  tested.  First,  we  model  telese- 
ismic  body-waves  of  two  events,  which  are  used  as  masters.  Short  period  depth 
phases,  pP  and  sP ,  are  used  to  establish  the  epicentral  depth,  and  the  events  are 
relocated  using  calibrated  stations  and  a  mantle  model  derived  for  this  region, 
TIP.  The  events  moved  upward  by  12  and  31  kms,  respectively.  The  regional 
waveforms  recorded  at  the  IRIS  station  GAR  from  the  best  determined  source  are 
forward  modeled  to  establish  a  local  crustal  model.  A  four-layered  model  with  a 
thickness  of  65  km  proves  effective.  Synthetics  from  this  crustal  model  can  then 
be  compared  to  the  data  from  other  events  where  the  depths  are  poorly  known,  i. 
e.,  PDE  depths  of  33  km.  Forward  modeling  of  the  short  period  phases  allows 
better  estimates  of  depth  and  makes  relocation  possible.  Next,  we  perform  a  long 
period  inversion  of  whole  and  partial  waveform  data  to  obtain  source  mechanism 
and  moment,  and  repeat  the  procedure  if  necessary.  This  technique  is  applied  to 
events  in  the  tectonically  active  Pamir-Hindu  Kush  region  to  test  its  usefulness. 
Five  crustal  events  sampling  various  azimuths  are  presented  as  examples  of  reloca¬ 
tion  and  determination  of  source  mechanisms  from  small  events.  The  smallest 
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event  has  a  moment  of  5.1  X  1022  dyne-cm.  For  the  seven  events  we  studied,  the 
average  depth  correction  to  PDE  is  19  km,  and  the  location  correction  is  11  km  on 
average.  The  method  can  be  used  to  identify  earthquakes  and  thus  it  lowers  the 
threshold  below  present  m4  :MS  discrimination,  since  these  events  are  too  small  to 
be  seen  teleseismically.  Secondly,  these  events  can  be  used  as  "masters"  in  the  cali¬ 
bration  of  other  systems,  case- based  event  characterizations,  etc. 

Running  title:  Source  Retrieval 

INTRODUCTION 

Earthquake  source  mechanisms  remains  an  important  topic  in  seismology. 
Techniques  to  obtain  accurate  seismic  parameters  from  teleseismic  waveform  data 
have  become  well  developed  both  for  surface  waves  (  Dziewonski  et  al.,  1981  )  and 
for  body  waves  (  Langston  and  Helmberger,  1975  ).  PDE  gives  the  source  mechan¬ 
isms  for  events  with  magnitude  greater  than  5.5  routinely  and  mechanisms  are 
published  in  bulletins  and  catalogs.  However,  smaller  earthquakes  in  the  magni¬ 
tude  range  4  <  mb  <  5.5  are  often  not  well  recorded  teleseismically.  The  impor¬ 
tance  of  source  mechanisms  of  these  smaller  earthquakes  is  growing  since  they 
prove  useful  in  tectonic  studies,  geodynamics,  earthquake  predictions,  and 
verification  of  test  ban  treaties.  However,  source  mechanism  for  these  small  events 
are  difficult  to  determine,  since  they  are  poorly  recorded.  Usually,  they  appear  on 
only  a  few  stations,  maybe  only  one  at  regional  distances.  In  this  paper,  we  ad¬ 
dress  the  issue  of  source  estimates  using  one  broadband  station  (  GAR  )  situated 
in  the  remote  region  of  Garm,  CIS,  see  Figure  1. 

GAR  is  one  of  the  IRIS  broadband  digital  stations  which  consist  of  the 
Widlandt-Streckeisen  sensor  and  the  Quanterra  data  logger.  The  systems  remain 
stable  for  even  the  rather  strong  motions  received  from  large  regional  events,  see 
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Figure  2.  This  feature  is  essential  to  our  procedure  since  we  can  study  these 
records  with  the  source  characteristics  determined  by  teleseismic  means.  Included 
in  Figure  2  are  the  responses  of  conventional  instruments.  For  example,  the  mo¬ 
tions  of  the  LP3090  are  similar  to  the  long-period  WWSSN  and  would  be  off-scale 
at  the  normal  gains  of  1500.  Source  inversion  studies  using  regional  phases  are 
usually  conducted  on  the  longer  periods  using  the  LP3090  as  a  filter  (  see  Dreger 
and  Helmberger  1991a,  1991b;  Fan  and  Wallace  1991  ).  The  short  period  Wood- 
Anderson  (  WA.SP  )  seismograms  usually  are  used  in  estimating  the  local  magni¬ 
tude,  Mi ,  and  prove  useful  in  defining  depth  phases.  A  strategy  for  determining 
a  crustal  model  that  produces  Green’s  function  explaining  regional  recordings, 
such  as  in  Figure  2,  is  given  in  a  earlier  paper,  Zhao  and  Helmberger  (  1991a  ). 
We  apply  the  same  procedure  in  this  paper  where  we  forward  model  broadband 
regional  waveforms  of  these  master  events  to  determine  the  velocity  structure.  We 
use  this  structure  to  predict  the  waveforms  of  the  smaller  events  where  we  gen¬ 
erally  find  disagreements  with  the  observations  in  ranges  and  depth  relative  to 
those  given  in  bulletins.  Once  the  source  depth  and  source  receiver  distance  of  a 
particular  earthquake  are  adjusted  to  fit  the  Green’s  functions,  we  use  the  inver¬ 
sion  code  of  Dreger  and  Helmberger  (  1991a,  1991b  )  to  determine  the  source 
mechanism.  The  physical  information  that  makes  this  method  work  is  extracted 
from  the  relative  P,  SV,  and  SH  strengths  as  pointed  out  by  Langston  (  1982  ). 

VELOCITY  STRUCTURE  BENEATH  THE  GARM  REGION 

In  this  section,  we  discuss  the  development  of  the  model  using  the  master 
events.  We  begin  by  first  determining  the  source  mechanisms  of  the  two  large 
events  displayed  in  Figure  2  and  their  detailed  source  properties.  Actually,  Event 
89205  will  be  used  as  the  master  event  and  Event  90036  will  be  used  in  a  compati¬ 
bility  test.  This  allows  an  independent  check  on  the  regional  results  compared 
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with  the  teleseismic  data,  and  provides  some  confidence  in  the  procedure. 

Locations  and  source  mechanisms  of  the  master  events 

In  order  to  get  more  accurate  information  about  these  sources,  we  forward 
model  the  teleseismic  waveforms  of  long-  and  short-period  waveforms  available 
from  GDSN  stations  to  determine  the  source  depth  and  source  mechanism.  Then 
we  assume  the  source  depth  to  be  fixed  and  use  the  relocation  scheme  of  Zhao  and 
Helmberger  (  1991b  )  to  determine  the  location  and  origin  time.  The  model  used 
to  calculate  Green’s  functions  is  TIP  (  Tibet  model,  Zhao  et  al.,  1991  ).  For  the  re¬ 
location  scheme,  a  TIP  source  region  and  Jeffreys-Bullen  receiver  regions  are  as¬ 
sumed.  Dziewonski  and  Anderson’s  (  1983  )  station  corrections  are  used  as  in  the 
above  study. 

Figures  3  gives  the  results  of  forward  modeling  the  89205  event.  The  event 
information  is  given  in  Table  1,  according  to  the  various  agencies.  In  Figure  3a, 
adequate  fits  for  the  long-period  waveforms  are  reached  for  all  stations  except 
WMQ.  WMQ  is  at  the  distance  is  14.9°,  within  the  upper  mantle  triplication,  and 
the  poor  fit  means  that  TIP  model  does  not  work  well  for  this  path.  TIP  works 
well  for  the  paths  crossing  the  Tibetan  Plateau  (  see  the  fits  of  KMI  and  LZH 
data,  Figure  3a  ).  The  short-period  synthetics  (  Figure  3b  )  indicate  the  triplica¬ 
tions  where  the  additional  arrivals  are  associated  with  pP  and  sP.  These  depth 
phases  are  seen  in  the  teleseismic  stations  ANTO,  COL,  CTAO,  HIA,  MDJ,  SSE 
and  TOL.  Our  forward  modeling  results  indicate  the  source  depth  is  85  km  (  PDE 
97  km  );  source  mechanism:  strike  245°,  dip  45°  and  rake  140°,  with  a  moment  of 
1.0  X  1025  dyne-cm.  The  location  is  36.131°N,  71.073°E  ,  and  the  origin  time  is 
3:27:48.9.  These  results  can  be  compared  with  the  other  estimates  in  Table  1. 

Forward  modeling  of  the  source  mechanism  of  the  90036  event  is  given  in 
Figure  4.  The  long-period  fits  are  not  as  good  as  those  for  the  event  89205  (  Fig- 
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ure  3a  ).  However,  good  fits  are  obtained  for  stations  CTAO,  GRFO,  KMI, 
MAJO,  and  WMQ  (  Figure  4a  ).  The  fits  are  also  acceptable  for  stations  ANTO, 
BCAO,  HIA,  NWAO,  and  SLR,  if  we  disregard  the  amplitude  ratio  of  P  and  pP 
waves.  There  are  also  three  arrivals  on  the  synthetics  (  except  the  stations  WMQ 
and  LZH,  Figure  4b  ).  Good  separations  of  the  three  arrivals  are  obtained  for  sta¬ 
tions  ARU,  BJI,  CTAO,  KMI,  LZH,  MAJO,  SLR  and  TOL,  and  as  well  as  the 
spacing  of  the  six  arrivals  at  stations  LZH  and  WMQ  which  fixes  the  source 
depth.  Results  from  forward  modeling  of  the  00036  event  indicate  a  source  depth 
of  100  km  (  PDE  131  km  ),  source  mechanism:  strike  115°,  dip  52°  and  rake  126°; 
with  a  moment  of  1.5  X  1025  dyne-cm.  The  relocation  yields  the  coordinates 
37.012*N,  7l.240°E,  and  the  new  origin  time  5:16:44.5.  These  results  can  be  com¬ 
pared  with  the  parameters  given  by  different  agencies  in  Table  1.  The  relocated 
events  agree  well  with  the  ISC  estimates  in  depth  and  origin  times  but  differ  con¬ 
siderably  from  HRVD  results. 

Velocity  structure  of  Garm  region 

Since  the  master  event  occurred  at  source  depths  not  usually  modeled,  we 
will  begin  with  some  preliminary  results  for  a  simple  layer  over  a  half-space 
model,  see  Figure  5.  The  responses  of  the  three  fundamental  fault  orientations  are 
displayed  to  gain  insight  into  possible  solutions.  The  SH  system  is  relatively  sim¬ 
ple  and  has  been  neglected.  The  upper  panel  shows  the  first  few  sets  of  ray  paths 
used  in  the  construction  of  the  synthetic  response  along  with  an  approximate  solu¬ 
tion  including  only  eight  rays,  the  direct  with  conversions  at  the  Moho  and  inter¬ 
nal  reflections  without  conversion  at  the  Moho.  With  only  eight  rays  to  consider, 
it  proves  relatively  easy  to  investigate  models  that  might  match  the  89205  obser¬ 
vations.  Figure  6  displays  the  best  fitting  one-layer  and  four  layered  models  found 
by  the  trial-and-error  method.  This  region  has  a  crust  of  about  65  km  thick,  and 
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a  mantle  with  compressional  velocity  of  about  8.0  km/sec  and  a  shear  velocity  of 
about  4.60  km/sec,  Table  2.  These  models  are  nearly  the  same  as  TIP  (  Zhao  et 
ah,  1991  )  and  thus  TIP  can  be  used  in  relocating  events  in  this  region. 

Synthetics  for  these  models  are  given  in  Figure  7.  Despite  the  complexity  of 
the  waveforms,  good  Gts  are  obtained  for  the  major  arrivals  of  P,  SP,  PPmP, 
SSmS  and  S  phases.  After  the  S  arrival  the  waveform  becomes  difficult  to  ex¬ 
plain  in  terms  of  a  flat-layered  model  and  will  be  neglected.  Note  that  the  arrivals 
SPmP  and  SP0  P  (  we  use  "o"  here  after  to  denote  internal  crustal  reflections  ) 
are  arriving  at  almost  the  same  time  on  the  data,  but  not  in  the  synthetics.  We 
could  find  a  one  layer  model  to  meet  the  requirement  of  all  the  internal  and  Moho 
reflections  arriving  at  about  the  same  time,  but  the  synthetic  fit  of  the  latter  part 
of  the  record  is  not  as  good  as  that  of  the  four  layer  model.  This  suggests  a  two 
dimensional  property  along  the  path. 

We  invert  the  whole  waveform  and  obtain  strike  272®,  dip  44°  and  rake 
136°,  a  mechanism  quit*  close  to  that  of  teleseismic  modeling.  Figure  8  gives  the 
comparison  of  long-period  WWSSN  data  and  synthetics  of  the  two  source 
mechanisms.  We  can  see  that  these  two  mechanisms  work  equally  well  in  match¬ 
ing  the  data.  If  we  do  not  know  the  source  mechanism,  we  have  no  reason  to 
choose  one  over  the  other.  The  moment  obtained  from  the  inversion  is  6.4  X  1024 
dyne-cm,  slightly  smaller  than  the  teleseismic  result. 

Relocation  by  forward-  modeling 

Most  of  the  smaller  events  are  not  well  constrained  in  depth,  as  indicated 
by  the  standard  PDE  default  of  33  km.  But  from  the  experience  gained  in  relocat¬ 
ing  events  beneath  Tibet,  Zhao  and  Helmberger  (  1991b  ),  we  expect  the  location 
tc  be  quite  good  in  comparison.  Essentially,  the  depth  trades  off  with  origin  time. 
This  means  that  we  must  explore  the  depth  and  to  a  lesser  extent  the  range  to 


47 


GAR  before  attempting  an  inversion.  Thus  we  compare  a  catalog  of  Green’s  func¬ 
tions  with  the  observed  waveforms  to  make  these  adjustments  (  Figure  9  ). 

Since  we  do  not  know  the  mechanism,  we  can  begin  by  exploring  with  the 
three  fundamental  faults  as  displayed  in  Figure  10.  The  three  dark  traces  are  the 
three  component  data  from  the  90064  event  and  the  eight  lighter  traces  are  the 
reflectivity  synthetics  for  the  fundamental  faults  (  three  traces  for  the  vertical, 
three  for  the  radial  and  two  for  the  tangential  ).  Since  any  fault  can  be  approxi¬ 
mated  by  a  linear  combination  of  these,  it  is  likely  that  one  of  these  traces  will 
look  like  the  data.  For  example,  the  synthetics  associated  with  the  45”  dip-slip  fit 
the  vertical  and  radial  components  of  data  quite  well.  For  the  tangential  com¬ 
ponent  the  synthetics  of  strike-slip  fault  seem  better,  which  is  expected  for  a 
thrust  mechanism.  We  compare  the  three  component  waveform  data  with  the  field 
of  synthetics,  by  assuming  the  source  mechanism  as  displayed  in  Figure  11,  and 
reach  a  distance  and  a  source  depth  that  produces  the  best  fitting  synthetics.  The 
distance  of  320  km  and  the  depth  of  10  km,  "320/10",  are  obtained  (  Table  1  ). 
Major  phases  are  lined  up  with  respect  to  PH  ,  the  depth  phase  sPmP,  Sn,  and 
surface  waves.  Next,  we  address  the  regioral  inversion  results  for  this  event  in 
conjunction  with  results  obtained  from  teleseismic  studies  since  this  is  a  relatively 
large  event. 

Event  90064 

Two  solutions  for  this  event  are  available,  namely,  strike  169®,  dip  30®, 
rake  —57®,  and  moment  2  X  1025  dyne-cm  from  NEIC;  and  strike  192®,  dip  36®, 
rake  —46®,  and  moment  2  X  1025  dyne-cm  given  by  Harvard,  see  Table  3.  The 
source  depth  is  3  -  60  km  from  different  agencies  (  see  ISC  bulletin  for  details  ). 

Figure  12  gives  the  comparison  of  data  with  the  synthetics  predicted  by  the 
different  source  mechanisms.  The  top  darker  traces  are  the  data.  The  second 


48 


traces  are  the  synthetics  of  the  source  mechanism,  strike  85°,  dip  77®,  rake  4®,  and 
the  moment  of  8.8  X  1024  dyne-cm  from  the  whole  waveform  inversion,  and  the 
third  trace  uses  only  a  portion  of  the  waveforms  denoted  by  arrows.  The  bottom 
traces  are  the  synthetics  of  the  source  mechanism  of  Harvard  CMT  solutions, 
strike  192®,  dip  36®,  rake  —  46®. 

Inversion  results  from  the  third  trace  are  very  close  to  the  CMT  solution 
and  are  almost  the  same  as  the  mechanism  given  by  NEIC.  This  solution  and  the 
CMT  solution  give  very  nice  surface  wave  fits  to  the  data  (  the  third  and  the  bot¬ 
tom  traces,  Figure  12  ).  However,  the  partial  waveform  solution  predicts  the 
wrong  polarity  of  both  Pn  and  Sn  (  SH  component  )  phases,  and  Harvard  CMT 
gives  a  nodal  Pn .  The  whole  waveform  solution  is  almost  a  pure  strike-slip,  how¬ 
ever  the  partial  waveform  solution  is  the  summation  of  pure  strike-slV,  ti.p-slip 
and  45®  with  about  same  strength  for  the  P—SV  waves.  The  source  mechanisms 
from  inversion  and  others  are  given  in  Table  3,  along  with  magnitude  informa¬ 
tion.  Thus,  from  Figure  12,  it  is  clear  that  the  source  mechanism  that  predicts  the 
best  fit  to  the  whole  waveforms  of  one  station  may  not  be  the  same  as  the  source 
mechanism  from  more  stations  and  points  out  the  limitation  of  the  one-station 
waveform  inversion.  We  will  address  these  issues  later. 

The  moment  obtained  form  inversions,  8.8  X  1024  from  the  whole 
waveform  data  and  1.1  X  1025  from  partial  waveform  data,  are  about  half  of  that 
given  by  NEIC  and  Harvard.  The  differences  may  be  caused  by  the  depth  uncer¬ 
tainty.  Note  that  the  amplitude  of  the  surface  waves  is  about  the  same  as  that  of 
the  5  waves  in  the  radial  component,  which  suggests  a  deeper  source,  but  this  is 
not  supported  by  the  vertical  component  of  the  data.  If  the  depth  is  10  km 
deeper,  one  can  expect  a  50  percent  increase  of  the  moment  (  Figure  9  ).  This  sug¬ 
gests  that  the  moment  is  an  uncertain  quantity  for  crustal  events,  if  the  depth  can 
not  be  established.  Next  we  return  to  the  other  master  event,  namely  the  deep 
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event  00036. 


Event  90036 

We  start  with  some  synthetics  generated  with  rays  containing  direct  and 
first  crustal  reflections  for  different  depths  and  distances  to  relocate  the  event.  The 
best  estimates  are  given  in  Figure  13  along  with  the  data.  Although  the  synthetics 
of  245/100  (  distance  245  and  depth  100  km  )  seems  better,  those  of  the  other  two 
are  also  acceptable.  The  depth  of  100  km  agrees  with  the  teleseismic,  for  the  verti¬ 
cal  travel  time  of  the  model  given  in  Table  2  is  almost  identical  to  the  model  TIP 
used  in  the  source  depth  determination.  The  distance  of  245  km  is  slightly 
different  from  that  of  teleseismic  results,  namely  235  km  and  that  of  PDE  230 
km.  However,  the  origin  time  of  the  relocation  is  1.7  seconds  earlier  than  that  of 
teleseismic  result  (  Table  1  ).  These  differences  can  be  easily  accommodated  by  an 
increase  in  mantle  velocity  beneath  85  kms. 

The  inversion  results  for  the  source  mechanism  for  the  event  90036  are 
given  in  the  lower  portion  of  Figure  13.  The  source  mechanism  is:  strike  227°,  dip 
47®  and  rake  46"  from  the  whole  waveform  inversion.  It  is  not  much  different  from 
that  of  the  teleseismic  result:  strike  245",  dip  50°  and  rake  53".  The  inversion 
results  are  displayed  as  the  second  traces  while  the  results  assuming  the  teleseism¬ 
ic  mechanisms  are  displayed  in  the  third  row.  The  mechanisms  are  quite  similar  as 
displayed  but  there  are  serious  problems  with  the  tangential  motions.  The  ob¬ 
served  tangential  motions  for  the  master  event  were  nodal  and  also  disturbed  in 
that  there  appears  to  be  (  P  -  5F  )  motions  observed  on  the  tangential  com¬ 
ponent.  No  clear  explanation  can  be  given  for  this  except  to  point  out  that  the 
strength  of  direct  S  is  greatly  reduced  by  the  transmission  across  the  mantle-crust 
boundary  at  these  angles,  that  is,  most  of  the  energy  is  reflected  down  into  the 
Earth.  Apparently,  the  non-planar  structure  becomes  particularly  important  in 
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this  geometry.  At  any  rate  the  lack  of  SH  control  greatly  reduces  the  effectiveness 
of  the  method.  Note  that  the  later  part  of  the  synthetics  is  about  the  same  for  the 
two  source  orientations,  whereas  the  teleseismic  solution  produces  a  slightly  small¬ 
er  P  -  wave  in  agreement  with  the  data. 

The  moment  from  the  inversion,  2.8  X  102S  dyne-cm,  almost  two  times  of 
that  of  teleseismic  studies,  1.5  X  1025,  is  the  same  as  that  given  by  NEIC  and 
Harvard  (  Table  3  ).  The  amplitude  predictions  of  the  teleseismic  mechanism  are 
smaller  than  that  of  the  data.  The  moment  would  be  1.8  X  1025,  instead  of 
1.5  X  1025,  in  order  for  the  teleseismic  source  mechanism  to  predict  the  correct 
amplitudes  of  this  regional  data.  Since  the  tangential  component  is  poorly 
matched,  we  did  not  include  it  in  the  moment  calculation. 

MORE  EXAMPLES 

Next,  we  present  more  examples  of  relocation  and  determination  of  source 
mechanisms  for  small  crustal  events,  89124,  89146,  90085,  and  91026.  Table  1 
gives  the  relocation  results  and  Table  3  gives  the  source  mechanisms. 


Event  89124 

The  forward  relocation  results  are  given  in  Figure  14,  in  which  the  event 
moved  from  PDE  208/33  to  200/20  .  Synthetics  at  small  perturbations  in  depth 
and  distance  are  included.  The  data  and  the  synthetics  are  lined  up  with  respect 
to  the  first  arrival  P0P  ■  If  we  use  the  depth  of  25  km  (  top  traces  ),  the  sPmP  ar¬ 
rival  is  later  than  in  the  data.  The  distance  of  210  km  predicts  a  separation  of  the 
P0P  and  SgS  phases  that  is  larger  than  in  the  data.  The  surface  wave  part  of  the 
data  is  very  complicated  on  the  vertical,  while  relatively  simple  on  the  radial  com¬ 
ponent.  The  SH  component  has  a  small  arrival  before  the  arrival  denoted  as 
SgS.  This  arrival  could  be  a  internal  reflection  in  the  crust,  and  suggests  that  this 


51 


particular  path  has  a  boundary  in  the  upper  20  km  of  the  -  «n*.  This  is  the  indi¬ 
cation  of  the  complexity  of  the  structure. 

The  inversion  results  are  given  in  the  lower  portion  erf  Figure  14,  where  the 
strike  33°,  dip  21®  and  rake  10®  is  obtained.  The  moment  is  3.8  X  1023  dyne-cm. 
The  Rayleigh  waves  are  excluded  in  this  inversion  because  of  the  complex  vertical 
component. 

Event  89146 

From  forward  modeling,  we  obtain  445/35  as  the  best  location,  which  is 
close  to  439/26,  the  estimate  given  by  PDE.  Figure  15  gives  the  synthetic  com- 
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parison  of  the  data  with  different  source  depths  and  distances.  The  data  and  syn¬ 
thetics  are  lined  up  with  respect  to  the  Pn  arrival  on  the  vertical  and  radial  com¬ 
ponents,  and  Sn  on  the  tangential.  For  the  vertical  and  radial  components,  the 
445/35  synthetics  fit  better  than  others.  However,  for  the  tangential  component, 
the  445/30  synthetics  are  better  (  better  timing  of  SmSSmS  ).  The  synthetics  of 
the  surface  waves  are  acceptable.  The  strong  Sm  SSm  S  phase  on  the  radial  com¬ 
ponent  suggests  a  relatively  sharp  Moho.  The  best  source  mechanism  obtained 
from  waveform  inversion  is:  strike  128®,  dip  89®  and  rake  —1®  and  the  moment  is 
1.4  X  1023  dyne-cm.  This  is  essentially  a  strike-slip  event.  In  the  lower  portion  of 
Figure  15  gives  the  synthetic  comparison  of  the  data  for  both  long-period  and 
broadband.  In  this  waveform  inversion,  we  used  the  445/35  synthetics  for  the 
vertical  and  radial  components  of  Green’s  function  and  the  445/30  synthetics  for 
the  tangential  component.  The  source  mechanism  from  partial  waveform  (  after 
Sn  arrival  )  inversion  is  the  same  as  the  results  given  above. 

Event  90085 

The  relocation  results  are  given  in  Figure  16  where  we  obtain  (  315/10  ) 
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compared  to  the  PDE  (  341/33  ).  Note  that  a  deeper  source  predicts  a  late  sPmP 
phase  (  top  traces  ),  larger  distance  gives  a  late  arriving  Airy  phase  (  bottom 
traces  ).  The  data  and  synthetics  are  lined  up  with  respect  to  Pn  phase  as  in  other 
plots.  Note  that  many  sharp  arrivals  on  the  synthetics  are  the  result  of  the  sharp 
boundaries  in  the  crust  of  the  model.  The  fundamental  faults  used  in  this  figure 
are  45®  dip-slip  for  the  vertical  and  radial  components  and  strike-slip  for  tangen¬ 
tial  component. 

The  source  mechanism  is:  strike  320®,  dip  24®,  rake  —  81®  and  the  moment 
5.1  X  1022  dyne-cm,  from  whole  waveform  inversion,  see  Figure  16.  This  source 
mechanism  is  very  similar  to  the  partial  waveform  result  of  Event  90064.  Howev¬ 
er,  a  clear  downward  polarity  of  Pn  shows  up  in  the  data.  This  event  is  small, 
and  high-level  long-period  noise  can  be  seen  on  both  radial  and  tangential  com¬ 
ponent.  The  source  mechanism  from  partial  waveform  (  after  Sn  arrival  )  inver¬ 
sion  is:  strike  316®,  dip  24®  and  rake—  83®.  It  is  almost  the  same  as  the  whole 
waveform  inversion. 

Event  91026 

The  relocation  results  for  this  event  are  given  in  Figure  17,  where  we  ob¬ 
tain  225/10  compared  with  the  PDE  225/32  estimate.  Four  major  arrivals  are  in¬ 
dicated  on  the  figure.  They  are:  P0P,  a  internal  crustal  reflection,  PmP ,  a 
reflection  from  Moho,  and  corresponding  critical  angle  reflections  of  5  waves.  The 
data  and  synthetics  are  lined  up  with  respect  to  the  first  arrival.  Slight  shifts  in 
distance  and  depth  indicates  that  the  225/10  is  the  best  model  as  indicated  in  the 
upper  portion  of  Figure  17. 

We  obtain  two  different  source  mechanisms,  namely,  strike  157®,  dip  83®, 
rake  95®,  and  moment  7.2  X  1023  dyne-cm  from  whole  waveform  inversion;  and 
strike  176®,  dip  16,  rake  288®,  and  moment  4.5  X  1023  dyne-cm,  from  a  partial 
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waveform  inversion.  Although  the  values  of  dip  and  rake  of  the  source  mechan¬ 
isms  are  very  different,  the  radiation  patterns  of  them  are  not  much  different  for 
this  particular  station.  The  synthetic  comparison  of  the  two  source  mechanisms 
are  given  in  the  lower  portion  of  Figure  17.  Amazing  fits  of  surface  waveforms  are 
reached  from  the  partial  waveform  source  mechanism,  which  predicts  the  wrong 
polarity  for  the  early  portion  of  the  data  (  bottom  traces  ).  The  whole  waveform 
solution  gives  a  good  fit  for  the  early  portion  of  the  data,  but  its  surface  wave  is 
not  as  good  as  that  of  the  partial  waveform  solution.  Comparing  the  waveforms 
of  the  data  and  synthetics  of  the  two  source  mechanisms,  we  conclude  that  the 
91026  event  has  a  nearly  dip-slip  source  from  the  vertical  and  radial  components. 
However,  a  strike-slip  source  predicts  a  better  fit  to  the  tangential  component. 

DISCUSSIONS 

Our  results  are  preliminary  in  nature  in  that  we  do  not  have  much  experi¬ 
ence  in  broadband  modeling,  given  the  limited  amount  of  recording  history. 
Nevertheless,  it  appears  that  modeling  regional  data  is  easier  than  generally  en¬ 
visioned.  In  fact,  some  of  our  fits  look  quite  promising.  In  this  section  we  will  dis¬ 
cuss  some  of  the  questions  of  uniqueness  and  possible  improvements  in  our  present 
method. 

The  Inversion  scheme  employed  in  the  above  examples  is  not  particularly 
sophisticated.  For  instance,  instead  of  applying  a  method  that  depends  on  the 
starting  model  we  could  have  employed  the  moment  tensor  approach,  see  Fan  and 
Wallace  (  1991  ).  There  is  also,  the  measure  of  fit  question,  Ll  versus  L2,  etc.  But 
since  we  are  dealing  with  very  approximate  Green’s  functions,  we  must  be  con¬ 
cerned  with  the  actual  fits  of  synthetics  to  observations,  and  a  number  of  parame¬ 
ters  involved.  Sometimes  the  surface  waves  match  very  well  but  S„  has  opposite 
polarity.  Since  the  phase  of  the  surface  waves  can  be  distorted  easily  by  shallow 
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structure,  we  may  want  to  choose  a  model  that  fits  the  Sn  polarity  even  though 
that  model  does  not  produce  the  best  inversion  criterion.  Thus  the  question  of 
selective  weighting  appears  where  body  waves  are  given  preference.  But,  perhaps, 
Pn  or  5n  is  nodal,  and  this  waveform  may  not  be  so  meaningful,  and,  perhaps, 
the  Love  wave  should  be  emphasized.  We  can  address  these  issues  to  some  extent 
by  conducting  numerical  experiments  in  which  we  examine  the  complete  source 
model  space  for  small  changes  in  the  Green’s  functions.  To  address  these  issues  of 
uniqueness,  we  choose  the  event  90064  as  a  test  case  since  it  has  a  CMT  solution 
available  for  comparison.  This  event  is,  also,  interesting  in  that  the  P- waves  are 
nodal  which  causes  difficulties  with  the  inversion,  as  discussed  earlier.  The  inver¬ 
sion  scheme  minimizes  the  error  defined  by 

T 

«,-/  I  /,(<)- 

0 

where  /,-,  g{  are  data  and  synthetics,  T  is  the  time  segment  used  in  the  inversion. 
Note  that  this  measure  of  fit  tends  to  emphasize  the  longer  periods  and  the  largest 
amplitudes.  In  Figures  18  and  19,  we  explore  the  possible  solutions  by  construct¬ 
ing  three  cross-sections  through  the  3D  source  space  in  terms  of  strike:rake, 
dip:rake,  and  strike:dip.  Synthetics  corresponding  to  the  best  fitting  models  are 
displayed  in  Figure  19.  The  top  panel,  Figure  18,  displays  the  parameter  space 
sampled  at  five  degree  intervals,  with  a  moment  of  10  X  1024  dyne-cm,  and  con¬ 
tours  in  dB  units  with  a  strike  of  350°,  a  dip  of  85°  and  a  rake  of  180°  yielding 
the  lowest  error.  The  synthetics  for  this  strike-slip  orientations  are  displayed  along 
the  traces  second  from  the  bottom  in  Figure  19.  The  other  solutions  discussed 
earlier  are  included  for  comparison.  The  normalized  cross-correlation  coefficients 
are  given  to  indicate  the  goodness  of  fit.  Note  that  e,  depends  on  the  amplitudes 
as  well  so  that  large  cross-correlation  coefficients  do  not  completely  control  the 
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solution  with  respect  to  the  inversion.  The  second  panel  in  Figure  18  displays  the 
parameter  space  generated  with  a  smaller  moment,  M0  =  5  X  1024.  The 
corresponding  synthetics  are  given  at  the  bottom  of  Figure  19  with  similar  results 
to  the  top  panel;  dropping  the  moment  tends  to  reduce  the  error  somewhat.  The 
third  panel  from  the  top  of  Figure  18  (  325/10  )  displays  even  smaller  when 
the  range  is  changed  by  5  kms.  The  rake  changes  substantially  and  the  overall  fits 
are  improved  both  in  cross-correlation  coefficients  and  reduced  overall  errors,  see 
the  top  set  of  traces  in  Figure  19.  The  bottom  panel  displayed  in  Figure  18  gives 
the  results  for  a  5  km  change  in  depth.  This  set  of  Green’s  functions  allows  the 
best  overall  fit  as  indicated  by  the  smallest  error  and  greatest  coefficients,  see  the 
second  set  of  traces  in  Figure  19. 

So  far  we  have  not  used  the  broadband  features  pf  these  seismograms  but 
have  concentrated  on  the  longer  period  motions.  Including  shorter  periods  or  actu¬ 
ally  modeling  both  short  period  and  long  period  simultaneously  is  quite  difficult. 
Recent  studies  of  TERRAscope  data  suggests  that  the  relative  timing  between 
phases,  i.  e.,  PmP,  SmS  etc,  while  not  so  noticeable  at  long  periods,  becomes 
obvious  in  short  period  seismograms.  Furthermore,  lateral  variation  in  structure 
causes  surface  reflected  depth  phases  to  show  ccLdiderable  scatter  which  can  be 
accounted  for  some  degree  by  applying  decomposed  Green’s  functions  in  which 
timing  shifts  between  phases  are  allowed  and  determined  by  local  calibration 
events  (  see  Helmberger  et  al.,  1992  ). 

A  less  sophisticated  approximation  that  still  uses  the  energy  ratios  of  short 
period,  Wood-Anderson,  over  long  period,  Press-Ewing,  band  pass  of  the  data. 
This  ratio  distinguishes  earthquakes  from  explosions  as  pointed  out  recently  by 
Woods  et  al.  (  1992  ).  The  ratio  appears  to  be  relatively  independent  of  distance 
and  depth  if  we  eliminate  the  surface  waves.  Adding  the  surface  waves  or  whole 
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record  ratios  show  a  dependence  on  depth  and  distance  (  Figure  20  ).  These 
features  can  be  seen  in  the  Green’s  functions,  Figure  9. 

Figure  20,  also,  displays  the  obvious  dependence  of  this  energy  ratio  on 
source  duration  or  stress-drop.  Thus,  we  can  obtain  a  rather  rough  estimate  of 
stress  level  by  adjusting  the  assuming  a  triangle  source  history,  to  match  the 
observed  ratio.  The  numbers  indicating  the  data  ratios  are  on  the  right  in  Figure 
19.  The  first  number  for  the  body  wave  portions  of  the  observed  and  synthetics 
agrees  roughly,  assuming  the  same  6t  =1.2  seconds  for  all  cases,  while  the  second 
number,  for  the  whole  seismograms,  is  effected  by  the  depth  and  mechanisms. 
Note  that  a  more  strike-slip  type  of  mechanism  reduces  the  the  long-period  sur¬ 
face  waves  and  increases  the  ratio,  bottom  traces.  A  deeper  source,  case  320/15, 
accomplishs  the  same  effect.  Thus,  including  such  constraints  can  help  establish 
the  source  depth  and  add  creditability.  However,  there  are  complications  which 
need  to  be  understood  and  individual  station  characteristics  established  before 
these  energy  ratios  can  be  treated  seriously.  For  example,  a  station  setting  on  soft 
material  will  probably  observe  anomalously  strong  short-periods.  Secondly,  direc¬ 
tivity  could,  also,  effect  these  energy  ratios.  These  issues  can  probably  be  assessed 
with  an  adequate  recording  history  and  local  experiences. 

Other  criterion  can  also  be  used  to  establish  the  nature  of  events  such  as 
the  relationship  of  the  derived  mechanisms  to  the  tectonics  of  the  region.  How¬ 
ever,  in  this  particular  stressed  zone  with  its  relatively  complex  geology,  see 
Chatelain  et  al.  (  1980  ),  we  probably  want  to  use  the  derived  mechanisms  to 
study  the  tectonics  as  opposed  to  deducing  mechanisms  from  the  knowledge  of 
particular  fault  dynamics. 

We  conclude  from  this  discussion  that  the  best  way  to  obtain  better  results 
is  to  obtain  more  accurate  Green’s  functions.  This  may  be  possible  by 
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regionalization  and  using  Green’s  functions  based  on  azimuth.  For  example,  in  the 
Pamir  region  (  latitudes  of  34  —  38°N  and  longitudes  68  —  72°E  ),  see  Figure  1, 
Roecker  (  1082  )  inverted  arrival  times  form  580  local  events  and  fixed  the  crustal 
thickness  at  70  km.  He  obtained  an  average  crustal  velocity  model  similar  to  ours. 
He  suggested  compressional  and  shear  velocities  of  8.08  and  4.63  km/sec  for  the 
depths  of  70-110  km.  These  results  are  in  good  agreement  with  the  model 
presented  in  Table  2  as  well  as  the  studies  to  the  south  by  Zhao  and  Xie  (  1992  ) 
and  Holt  and  Wallace  (  1990  ).  The  latter  two  studies  suggest  considerable  lateral 
variations  as  originally  proposed  by  Roecker  (  1982  ).  This  feature  is  well  observed 
in  a  deep  seismic  sounding  profile  running  from  (42®N,  72®E)  to  (34°N,  75°E), 
Beloussov  et  al.  (  1980  )  and  Kaila  (  1981  ),  where  the  crustal  thickness  changes 
from  50-55  km  at  the  northern  end  to  60-65  km  at  the  south  end.  This  means 
that  the  one  dimensional  velocity  model  assumed  in  our  analysis  is  a  rough 
approximation  and  allowing  some  regionalization  could  be  quite  helpful. 

CONCLUSIONS 

In  conclusion,  we  have  examined  the  possibility  of  providing  source  location 
and  characterization  using  a  single  modern  station.  We  selected  the  Hindu-Kush 
region  in  this  pilot  study  largely  because  of  it  high  seismic  activity,  numerous 
observations  with  only  a  short  recording  history.  Unfortunately,  a  high  rate  of 
seismic  activity,  especially  the  dip-slip  style  in  this  region,  is  associated  with  com¬ 
plex  crustal  structure  and  wave  propagation.  Nevertheless,  we  found  excellent 
waveform  fits  for  some  small  events  that  can  not  be  studied  teleseismically.  But 
considerable  liberty  was  taken  in  moving  the  sources  around  to  find  a  Green’s 
function  that  worked.  Such  adjustment  shift  the  timing  and  phases  strengths  of 
the  various  wave  packets.  Shifts  of  this  sort  could  be  caused  by  lateral  variations 
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in  structure  and  maybe  not  correspond  to  the  true  source  position.  However,  since  the 
applied  shifts  in  this  study  produced  solutions  comparable  to  the  teleseismic  solutions 
for  the  large  events,  we  think  the  results  should  also  apply  to  the  smaller  events.  To 
confirm  this  conclusion  requires  more  information,  perhaps,  some  short  period  regional 
data  (  networks  )  or  preferably  more  IRIS  stations. 
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